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EVALUATION OF ALTERNATIVE PROCEDURES FOR 
ATMOSPHERIC ABSORPTION ADJUSTMENTS 
DURING NOISE CERTIFICATION 


VOLUME IIT: TABLES OF ATMOSPHERIC-ABSORPTION LOSSES 


INTRODUCTION 


While analytical methods are available to calculate the attenuation of 
the amplitude of a freely propagating sound wave because of absorptive prec- 
esses in the atmosphere, for ready comparison of the magnitude of the dif- 
ferences in the results of applying the various methods it is more convenient 
to use tables that present data for a wide range of the relevant physical 
variables. The purpose of this volume is to present such tables for a range 
of air temperatures, relative humidities, frequencies, calculation methods, 


and pure-tone or broadband sounds. 


The quantity that is listed in the tables is the attenuation, in decibels, 


of the amplitude of a sound wave after propagation over a distance of 300 m. 


The attenuation values in the tables consider atmospheric-absorption losses 
only; attenuation caused by spreading the acoustic power over increasingly 
larger surface areas as the wave propagates away from the source is not 
included and would be additive to the tabulated attenuation values. A dis- 
tance of 300 m was chosen so that most of the tabulated low-frequency atten- 
uation values would be significant with a precision of 0.1 dB. The 300-m 
distance also closely corresponds with the 1000-ft (304.8-m) distance used 
with some published tabulations of pure-tone atmospheric-absorption coeffi- 


cient (i.e., absorption loss per unit distance) in decibels per 1000 feet. 


Two methods of calculating atmospheric absorption for pure-tone sounds 
are included: that of American National Standard ANS S1.26-1978 and that 
of the Society of Automotive Engineers Aerospace Recommended Practice 
ARP866A. Broadband sounds with band-level spectral slopes of +1, -6, and 
-12 dB/band are included for determination of attenuation by the band- 


integration method. The broadband sounds are considered to be analyzed 


by ideal 1/3-octave-band filters. The method of SAE ARP866A was used 
exactly as though it was being applied to a broadband sound analyzed by 
1/3-octave-band filters. The frequencies of the pure-tone sounds for which 
the attenuation was calculated by the method of ANS S1.26-1978 were those 
of the exact center or geometric mean frequencies of the ideal 1/3-octave- 


band filters. These choices are consistent with the analyses in Volume I. 


The remainder of this Volume consists of five sections. Sectidén 2 
presents a discussion of the analytical basis for the tables. Section 3 
offers some general remarks relative to possible uses for the tables. 
Section 4 describes the scope and content of the tables. Section 5 contains 
the tables themselves. Section 6 is a listing of the FORTRAN statements for 


the computer program that generated the tables. 
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2. ANALYTICAL BASIS FOR THE TABLES 


Attenuation values tabulated in Section 5 were calculated by five 
methods [designated as methods (A), (B), (C), (D), and (E)] for each 
combination of air temperature and relative humidity. For methods 
(A) to (D), atmospheric absorption was calculated using the analytical 
model in American National Standard ANS S1.26-1978. For method (E), 
the atmospheric absorption and the attenuation over the 300-m path were 
calculated by the method of SAE ARP866A. The analytical descriptions in 
this Section partially duplicate those in Volume I and are provided here 


for the reader's convenience. 
2.1 Method (A) 


Method (A) determines the attenuation due to atmospheric absorption 
over the 300-m path for pure-tone sounds at the exact center frequencies 
of the standard 1/3-octave-band filters. For a pure-tone sound, the sound 
pressure level, Lo, at a point 2 is related to the sound pressure level, 
Ly)» at point 1 by 


where (L, - L,) = 10 log P,*/P,” 
= -10 log P,*/p,* (1) 


with point 2 being farther from the source than point 1 and on a line 
in the direction of propagation. 


Neglecting losses caused by the geometric spreading of sound power 


over larger surface areas, the squared pressures are related by 
-20 
Py? 2 p,? en 208 (2) 


where a is the sound atmospheric absorption coefficient in nepers/m and 
€ 1s the propagation distance in meters. 


With eq. (2), the band level difference for a pure-tone sound is 


(L) - L,) = 20 a€ logig e 
= al (3) 


where a = (20 log e)(a) = 8.6859a is the atmospheric sound absorption coef- 
ficient in dB/m. 


For specified air temperature, relative humidity, air pressure, and 
frequency, the valueof a is calculated and multiplied by the propagation 
distance € = 300 m to determine the magnitude of the attenuation in decibels. 


2.2 Methods (B), (C), and (D) 


For methods (B), (C), and (D), the sound was assumed to be a broadband 
noise with a continuous and constant-slope spectrum. With G(f) representing 
the sound pressure spectral density function, the time-averaged mean squared 


pressures in a filter band at locations 1 and 2 are determined from 


and 9 
Po" =(" G(£) dé. (5) 
0 
At any frequency, the pressure spectral density functions are related 
by 
= -206 
G, Gg, e 
or, in an equivalent form, by 
2 -aE/10 
Gy (G,) (10 ) 
=G_A (6) 


1 
where A is the atmospheric absorption function for specified air temperature, 


relative humidity, and air pressure. 


By analogy with eq. (1), the band level difference is, with eq. (6) 


for Go» 
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(Ly L,) 10 log py /Py 
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Assuming that the filters have ideal transmission-response characteris~ 
tics, the infinite limits for the integrals in eq. (7) can be replaced by 
integration over the passband of the filter from the lower, fi to the 


upper, fy bandedge frequencies. 


For a sound pressure spectral density with a constant spectral slope, 


the noise spectrum at point 1 (the G, function) can be written as 
SL 
G, (£) K(£/f ) (8) 


where K is an arbitrary constant that can be taken outside the integrals, 
f. is the band center frequency, and SL is the slope. The constant K 
actually represents the pressure spectral density at the band center 
frequency fA its value is arbitrary since it appears identically in 


numerator and denominator. 


Using eq. (8) and the assumption of ideal filters, the band~level 
difference (or attenuation) becomes 
f £ 
a6 U SL U SL 
(Ly = Ly) = -10 roe|| f (£/£,) a atl/[/ (£/£,) ee|t. (9) 
f 


L L 
The denominator term in eq. (9) can be evaluated directly to give 


£ 
U SL 
f (£/£.) df = fo ln (£,,/£ 


D (10) 
fr 
when SL = -1, and 
f 
U SL SL +1 SL +1 
f° crey® at =, AGL + DITE,/ED - (£,/6,) } (1) 
fL 
when SL # ~1. 
Introducing the band frequency ratio 
RF fy/ fh (12) 
and noting that 
a 1/2 = 1/2 
f/f, RF and f/f, RF- (13) 


and restricting the slope SL to be not equal to -1, eq. (11) for the 
denominator term becomes 
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£ 
DEN =| u c£/£ df= £ [1/BS] [RFP9/2 — pe@BS/2y (14) 


*L 


where the notation BS has been used for the band slope (i.e., the slope of 


the 1/3-octave-band spectrum at point 1). The abbreviation DEN is used in 
the statements for the computer program. 


The band slope is related to the slope of the sound pressure spectral 
density by 


SL = BS - 1 


or SL + 1 = BS. 


The band slope BS is determined from the band-level differences using 
BS = AL/(10 log RF) (15) 


where AL is the band-level difference in dB/band. 


For 1/3-octave-band filters, the frequency ratio is 


RF = 10/19, (16) 
The band center frequency is found from 
f* 10 15BN/10 (17) 


where ISBN is the set of International Standard Band Numbers for 1/3-octave- 
band filters that ranges from 17 to 40 for nominal band center frequencies 
ranging from 50 to 10,000 Hz. 


With fo RF, SL (and hence BS) known or specified, the denominator 
term can be calculated readily using eq. (14). 


The numerator term in eq. (9) is evaluated by a numerical integration 
method with the limiting frequencies fh and fy defined by eqs. (12), (13), 


and (17) and with the attenuation function defined by eq. (6) for a propa- 
gation distance of € = 300.0 m. 
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a 


The attenuation, for each filter band, is then calculated by carrying 
out the operation defined by eq. 


+1, -6, and -12 dB/band. 


(9) for specified band-level slopes AL of 


Note that AL = +1 implies a white-noise spectrum 
(SL = 0) with constant pressure spectral density. 


A band-level slope of 
AL = 0 would mean a pink-noise spectrum with SL = -l. 
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3. REMARKS ON USES FOR THE TABLES 


The tables in Section 5 permit comparative analyses of attenuation 

values determined using the pure~-tone atmospheric-absorption calculation 
method of ANS $1.26-1978 at the exact 1/3-octave~band center frequencies. 

For broadband sounds with constant spectral slopes, attenuation values are 
calculated by a numerical integration scheme over the passband of corres- 
ponding ideal 1/3-octave-band filters with atmospheric absorption specified 

by $1.26-1978. Attenuation values calculated using the atmospheric-absorption 
model of SAE ARP866A are also provided. Tabulated values are given for a ' 


wide range of air temperature and relative humidity so that the magnitude 


of the differences among the calculation methods and various trends can be 


illustrated. 


Attenuation values were selected for the tables, and not atmospheric 
absorption coefficients, because the concept of an absorption coefficient 
is not directly applicable to a broadband sound analyzed by fractional- 
octave-band filters. Attenuation (the difference in band sound pressure 
levels at two points) depends on the length of the propagation path, the 
slope of the sound pressure spectral density at the first point on the 
path, as well as the atmospheric conditions along the path and the frequen- 
cies within the filter passband. Equations (6) and (7) in Section 2 show 
the essential linkage of these variables in a calculation of band-level 


attenuation. 


Attenuation values for method (E) were calculated by a straightforward 
application of the procedures of SAE ARP866A, including the shift from the 


use of nominal band center frequencies to nominal lower bandedge frequencies 


for the last four bands. This choice was made to keep the comparisons t 


consistent with current usage of ARP866A even though the resulting 
comparisons include a built-in distortion for band center frequencies of 
5000, 6300, 8000, and 10,000 Hz. 


The tables were prepared with the assumption that the air pressure 
was 1.0 standard atmospheres. Separate calculations have shown that air 
pressure is a minor variable over the range of air pressures usually 


encountered in aircraft noise measurements. 
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Good agreement between the attenuation calculated at the band-center 
frequency and that calculated over a filter passband by a band-integration 
method depends on the slope of the sound pressure spectral density not being 
too steep and on the fact that 1/3-octave-band filters are relatively narrow 
for low to mid frequencies. At high frequencies (or over long pathlengths) 
there can be relatively large differences between the attenuation at the 


band-center frequency and that calculated by a band-integration method. 


For high-frequency noise spectral slopes that are relatively steep 
(because of highly absorptive conditions, or long propagation paths, or both), 
inspection of the tables shows that the band-level adjustment factor from 
actual to acoustical reference conditions of 25°C and 70 percent relative 
humidity is generally smaller by the band-integration method than by either 
pure-tone method (A or E), at least for band center frequencies less than 
5000 Hz. This observation is consistent with the results reported in Volume 
I. 
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4. DESCRIPTION OF THE TABLES 


For each combination of air temperature and relative humidity, the 
300-m attenuation values in the tables are calculated by five methods for 
each frequency or frequency band. The methods are labeled (A), (B), (C), 
(D), and (E) over the columns in the tables. The various methods are 


described below: 


(A) The attenuation experienced by a pure-tone sound having a 
frequency equal to the exact center frequency of the 24 1/3-octave bands 
and with atmospheric absorption coefficients calculated by the method 
of American National Standard ANS S1.26-1978, 


(B) The attenuation, over the frequency range of the passband of 
ideal 1/3-octave-band filters, experienced by a broadband sound having a 
band-level slope of +1.0 dB/band, i.e., a white-noise signal. Atmospheric 
absorption calculated by ANS S1.26-1978 at specified frequencies over the 
range of each filter passband. Attenuation is calculated by a numerical 


integration process. 


(C) Like (B), except for a band-level slope of -6 dB/band. 


(D) Like (B) and (C), except for a band-level slope of -12 dB/band. 


(E) The attenuation as calculated by the method of the Society of 
Automotive Engineers’ Aerospace Recommended Practice, ARP866A, 15 March 
1975, for broadband noise signals analyzed by 1/3-octave-band filters. 
Note that the method is really a calculation of the atmospheric-absorption 
loss of a pure-tone sound but, in contrast to method (A), the calculation 
is done for the nominal band center frequencies for bands 1 to 20 and for 
the nominal lower bandedge frequencies for bands 21 to 24 (i.e., nominal 
band center frequencies from 50 to 4000 Hz and nominal lower bandedge 


frequencies for band center frequencies from 5000 to 10,000 Hz). 
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| 5. TABLES OF ATMOSPHERIC ABSORPTION LOSSES 


The tables on the next 34 pages present values of the attenuation, 
in decibels, caused by atmospheric absorption losses for a sound wave 
propagated over a distance of 300 meters. Each page is for a specified 
air temperature ranging from 2 to 35 degrees celsius at intervals of one 
degree. Each page shows attenuation values for relative humidities ranging 
from 10 to 100 percent at intervals of 10 percentage points. For each 


combination of air temperature and relative humidity, attenuation values 


a are presented for 24 frequencies, or frequency bands, corresponding to 
iz the 1/3-octave bands having nominal band center frequencies ranging from 
b | 50 to 10,000 hertz. The air pressure is assumed to be 1.0 standard 


atmosphere for all calculations. 


Column headings (A), (B), (C), (D), and (E) for the methods of calcu- 
lating attenuation are defined in Section 4. 
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6. 


LISTING OF PROGRAM STATEMENTS 


The next 12 pages list the statements in FORTRAN IV for the computer 


follow the analytical development in Section 2 for methods (A) to (D). Calcu- 
lation of attenuation by SAE ARP866A in method (E) essentially uses subroutine 
ARP866 from Volume II. Numerical integration over the filter passband for 
| methods (B), (C), and (D) uses the standard method in subroutine QSF from 
the IBM Scientific Subroutine Package. This is the same numerical integration 


| 
‘ | program that generated the tables in Section 5. The computation algorithms 


{ 
| subroutine used in Volume II; it is included here to provide a complete 
1 


package. 


The calculation process follows the format selected for the tables and 


i is accomplished in two steps for each air temperature. In the first step, 


attenuations by the five methods are calculated for five relative humidities 


from 10 to 50 percent and for 24 frequency bands. The second step repeats 


the process for the second set of five relative humidities from 60 to 100 


percent. The entire process is then repeated for the next increment in air 


temperature until the calculations have been completed for all 34 tempera- 


tures. 


The program is self contained; no input data are required for execution. 


A number of descriptive comments have been included to facilitate the incor- 


poration of modifications, if desired, and to clarify the various steps in r 


the calculations. 


1.000 C##% PROGRAM NAME IS TABLES ' 


| 2.000 C 
7 3.000 Re 38 38 36 98 36 38 38 SE 36 56 SE SE SE SE OE SE SE SE SE SE SE SE 38 OE SC SE SE SE 98 86 36 
4.000 C 
| 5.000 C PROGRAM WAS PREPARED BY DYTEC ENGINEERING, INC., OF 
| 6.000 C LONG BEACH, CA; VERSION 2 JUNE 1979. 
7.000 C 
| 8.000 C PURPOSE OF THE PROGRAM IS TO CALCULATE AND LIST TABLES FOR THE 
9.000 C ATTENUATION OF SOUND OVER A PROPAGATION DISTANCE OF 300 METERS. 
10.000 ¢ 
11.000 Cc ABSORPTION COEFFICIENTS ARE CALCULATED BY THE METHODS OF 
12.000 C AMERICAN NATIONAL STANDARD ANS $1.26-1978 AND BY 
13.000 C SAE AEROSPACE RECOMMENDED PRACTICE ARP866A. 
14.000 C 
15.000 C ATTENUATION VALUES ARE CALCULATED FOR 34 AIR TEMPERATURES, 
16.000 C (2 TO 35 DEG C AT 1 DEG C INTERVALS), 10 RELATIVE HUMIDITIES (10 TO 
17.000 C 100 PERCENT AT 10 PERCENTAGE-POINT INTERVALS), AND 24 FREQUENCIES 
; 18.000 C CCORRESPONDING TO THE CENTER FREQUENCIES OF THE 1/3-OCTAVE BANDS 
19.000 C BETWEEN 50 AND 10,000 HZ). AIR PRESSURE IS ASSUMED TO BE CONSTANT 
20.000 C AT 1.0 STANDARD ATMOSPHERE. 
21.000 ¢ 
22.000 C THE PROGRAM IS SELF CONTAINED. NO INPUT DATA ARE REQUIRED. 
23.000 C 
24.000 C OUTPUT IS ON DEVICE 6 WITH 58 LINES PER PAGE. 
25.000 C 
26.000 C ATTENUATION VALUES ARE LISTED IN DECIBELS OVER A 300-M 
27.000 C OISTANCE. 
28.000 C 
29.000 C28 FE 8 08 a8 8 28 BE 8 28 ok BE 8 08 BE SE OE OE 8 28 OE 8 38 OS 38 28 28 BE 28 OE 28 ES 28 28 58 28 OE 8 OE 28 26 28 28 28 8 28 OS 28 8 oe 2 SE oS eae eh ce sega ae cece ce 
30.000 C 
31.000 DIMENSION FREQ(24), IFREQC24),$126(10), 
32.000 1SL1P(10),SL6NC10),SL12N(10), 
33.000 2SAE866(10),STEPC24), TABLE(58),Y1(31), 
34.000 3¥2€31),Y3C31),Z1031),22¢€31),23¢€31) 
35.000 C 
36.000 LOGICAL FRQTIC51) 
37.000 C 
38.000 REAL ISBN,LPSOPO 
39.000 C 
40.000 DATA IFREQ /50,63,80,100,125,160, 200, 250, 315, 400, 500, 630, 800, 
41.000 11000, 1250, 1600, 2000, 2500, 3150, 4000, 5000, 6300, 8000, 10000/ 
42.000 C 
43.000 DATA FRQTI/' ',!' ',!' #,! PONE tOtG TM TINS TAT IES, 
44,000 1'O','N','E!,! PTS VAS OPER 3 tO! MBS TAY NGS TD, 
45.000 2! bo TES NG, fT, tet Re TS OR! Se Oe US SEN, 
' 46.000 SCE YE NS ht tA pt Z ie 4 '] 
47.000 C 
48.000 DATA TABLE /0.00,0.00,0.250,0.315,0.50,0.70,0.6,0.84,0.7,0.93, 
49.000 10.8,0.975,0.9,0.996,1.0,1.0,1.1,0.97,1.2,0.9,1.3,0.84,1.5,0.75, 
50.000 21.7,0.67,2.0,0.57,2.3,0.495,2.5,0.45,2.8,0.4,3.0,0.37,3.30,0. 33, 
51.000 33.6,0.3,4.15,0.26,4.45,0.245,4.80,0.23,5.25,0.22,5.7,0.21, 
52.000 46.05,0.205,6.5,0.2,7.0,0.2,10.0,0.2/ 
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DATA FREQ /50.,63.,80.,100.,125.,160.,200.,250.,315.,400.,500., 
1630.,800.,1000.,1250.,1600.,2000.,2500.,3150.,4000.,4500., 
25600.,7100.,9000./ 


DATA STEP/13*8.,10.,12.,14.,16.,18.,20.,22.,24.,26.,28., 30./ 


FTE RT ERR OR ET ER wren 


DEFINE AIR PRESSURE, PA, IN STANDARD ATMOSPHERES 
PA=1.0 
DEFINE REFERENCE AIR TEMPERATURES 


T0=293.15 
T01=273.16 


DEFINE RF, THE EXACT VALUE OF THE FREQUENCY RATIO FOR 
1/3-OCTAVE BANDS. 


AAAN 


RF=10.0%*0.1 

FOR BAND-LEVEL DIFFERENCES OF +1.0 DB/BAND CI.E., A WHITE 
NOISE SPECTRUM), THE BAND-LEVEL SLOPE, SLOPE1, IS: 

SLOPE1=1.0/(€10.0*LOG10CRF)) 


AND THE CORRESPONDING SOUND PRESSURE SPECTRAL DENSITY SLOPE, 
SLI, IS: 


Cc 
Cc 
Cc 
Cc 
Cc 
c 
c 
Cc 
Cc 


SLI=SLOPE1-1.0 


FOR BAND-LEVEL DIFFERENCES OF -6.0 DB/BAND, THE BAND-LEVEL 
SLOPE, SLOPE2, IS: 


SLOPE 2=-6.0/(€10.0*LOG1O(RF)) 


AND THE CORRESPONDING SOUND PRESSURE SPECTRAL DENSITY 
SLOPE, SL2, IS: 


aAaNAAaAN ANAA 


SL2=SLOPE2-1.0 


FOR BAND-LEVEL DIFFERENCES OF -12.0 DB/BAND, THE SLOPES ARE: 


aAaAN 


SLOPE 3=-12.0/(€10.0*LOG10CRF)) 
SL3=SLOPE3-1.0 


INITIAL VALUES OF AIR TEMPERATURES IN DEG C AND DEG K 


TC=1 . 0 
TK=274.15 


Fe ENE RD PROT PRY LOO PTET ENET CPt 


105. 
106. 
107. 
108. 
109. 
110. 
111. 
112. 
113. 
114, 
115. 
116. 
117. 
118. 
119. 
120. 
121. 
122. 
123. 
124. 
125. 
126. 
127. 
128. 
129. 
130. 


131 


132. 
133. 
134. 
135. 
136. 
137. 
138. 
139. 
140. 
141. 


142 


143. 
144, 
145, 
146. 
147. 
148, 
149.000 
150.000 
151.000 
152.000 
153.000 
154.000 
155.000 
156.000 
157.000 
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000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 


aAaAAAN fo) a aAaANn a 


aAnAAN 


aAANAAN 


a anNAaAN 


10 


20 


30 


40 


50 


DO 270 [=1,34 


TC=TC+1.0 
TK=TK+1.0 


WRITEC6,10) TC 

FORMAT(8/,T28, ‘ATTENUATION, IN DB, OVER 300-M PATHLENGTH', 

1' FOR AN AIR TEMPERATURE OF ',12,' DEG C') 

WRITEC6, 20) 

FORMAT(T23,'AIR PRESSURE = 1.0 STD. ATMOSPHERE; RELATIVE’, . 
1' HUMIDITY, IN PERCENT, NOTED ABOVE COLUMNS') ; 
WRITEC6, 30) 

FORMATC/,1T17,'METHODS CA), (8), CC), CD), AND CE) FOR', 

1' CALCULATING ATTENUATION ARE EXPLAINED IN! 

2' ACCOMPANYING TEXT') 

WRITEC6, 40) 

FORMAT(/,T17,'10 PCT RH',16X,'20 PCT RH',16X,'30 PCT RH', 

116X,'40 PCT RH',16xX,'50 PCT RH') 

WRITEC6,50) 

FORMATCT8,5' CA) (CB) (CC) CD) CED') 


DO 140 J=1, 24 

INITIAL VALUE OF RELATIVE HUMIDITY, IN PERCENT 
HR=0.0 

DO 120 K=1,5 

HR=HR+10.0 


CALCULATE THE EXACT BAND CENTER FREQUENCIES FROM THE 
INTERNATIONAL STANDARD BAND NUMBERS, ISBN. 


ISBN=J+16 
FC=10.0**(C1SBN/10.0) 


THE EXACT LOWER, FL, AND UPPER, FU, BANDEDGE FREQUENCIES FOR 
THE FILTER BAND ARE: 


FL=FC/SQRTCRF) 
FU=FC#SQRTCRF) 


THE WIDTH OF THE EQUALLY SPACED FREQUENCY INTERVALS, DELF, 
OVER THE PASSBAND OF THE FILTER IS: 


DELF=CFU-FL)/STEPCU) 


THE VALUE OF THE FREQUENCY AT ONE STEP BELOW THE LOWER 
BANDEDGE FREQUENCY IS: 


F=FL-DELF 


ee re Carne 


158. 
159. 
160. 
161. 
162. 
163. 
164. 
165. 
166. 
167. 
168. 
169. 
170. 
171. 
172. 
173. 
174. 
175. 
176. 
177. 
178. 
179. 
180. 
181. 
182. 
183. 


184 


185. 
186. 
187. 
188. 
189. 
190. 
191. 
192. 
193. 
194, 
195. 
196. 
197. 
198. 
199. 
200. 
201. 
202. 
203. 
204, 
205. 
206. 
207. 
208. 
209. 
210. 


000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 
000 


AAANAANADAN AAD 


aan aan 


Aan 


ANAND 


AAAND aaa 


THE TOTAL NUMBER OF FREQUENCY STEPS IN THE NUMERICAL 
INTEGRATION OVER THE PASSBAND OF THE IDEAL FILTER IS: 


NN=L4+STEPCJ) 

CALCULATE THE DENOMINATOR TERMS USED IN THE INTEGRATION 
METHOD TO DETERMINE THE BAND ATTENUATION OVER EACH OF THE 
THREE SPECTRAL SLOPES. 

THE DENOMINATOR TERM FOR THE +1.0 DB/BAND SLOPE IS: 


DENL=FC#C1.0/CSL141.0))*CCRF#*CCSL141.0)/2.0))- 
1 CRF¥¥(-(§L14+1.0)/2.0))) 


THE DENOMINATOR TERM FOR THE -6.0 DB/BAND SLOPE IS: 


DEN2=FC#(1.0/CSL2+1 .0))*CCRF**CCSL 241 .0)/2.0))- 
1 CRFS*(-(CSL241.0)/2.0))) 


THE DENOMINATOR TERM FOR THE -12.0 DB/BAND SLOPE IS: 


DEN3=FC#(1.0/(SL341.0))*CCRF**CCSL341.0)/2.0))- 
1 CRF®*(-($L341.0)/2.0))) 


CA) CALCULATE ABSORPTION COEFFICIENTS BY $1. 26-1978 


LPSOPO0=10.79586*C(1.-CTO1/TK)) - 5.02808*LOG10CTK/TO1) 


1 +1.50474E-45(1.-10.%*C-8. 29692*%CCTK/T01)-1.))) 
2 +0.42873E-3%(10.**(4. 76955%C1.-CTO1/TK)))-1.) 
3 -2.2195983 


PSOP0=10.**LPSOPO 

H=HR=*=PSOP0/PA 
FRO2=PA#(24.44.,41E04*H*((0.05+H)/C0.3914+H))) 
FRN2=CPA/SQRTCTK/T0))*(9.+ 


1 350. *H#EXPC-6.142*CCCTK/T0)**(-1./3.))-1.))) 
ALPHA=CFC#*2.)*(C C1. 84E-11)*C1. /PA)*SQRTCTK/TO)) 

1 +CCTK/TO)#*(-5./2.))*#CC1. 278E-2% 

2 CEXPC-2239.1/TK))/CFRO2+CCFC##2.)/FRO2))) 

3 +€0. 1068 *CEXPC-3352./TK))/CFRN2+CCFC#*2.)/FRN2))))) 


THE ABSORPTION COEFFICIENT, AFC IN D8/M, AT THE EXACT SAND 
CENTER FREQUENCY BY S$1.26-1978 IS: 


AFC=8.686*ALPHA 
THE ATTENUATION OVER A 300-M DISTANCE IS: 
$126CK J=AFC#300.0 


CB), €C), AND (D0) CALCULATE BAND ATTENUATION FACTORS BY 
THE INTEGRATION METHOD FOR THE THREE SPECTRAL SLOPES. 
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211.000 
212.000 
213.000 
214.000 
215.000 
216.000 
217.000 
218.000 
219.000 
220.000 
221.000 
222.000 
223.000 
224.000 
225.000 
226.000 
227.000 
228.000 
229.000 
230.000 
231.000 
232.000 
233.000 
234.000 
235.000 
236.000 
237.000 
238.000 
239.000 
240.000 
241.000 
242.000 
243.000 
244.000 
245.000 
246.000 
247.000 
248.000 
249.000 
250.000 
251.000 
252.000 
253.000 
254.000 
255.000 
256.000 
257.000 
258.000 
259.000 
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ANAAND ANNAN ANnoOA 


ANAND ANAA 


60 


AANMAaAN 


aaa 


DO 60 JUJs1,NN 


DETERMINE THE FREQUENCIES, F, AT EACH STEP OVER THE 
FREQUENCY RANGE OF THE FILTER PASSBAND. 


F =F +DELF 


THE VALUES OF THE NORMALIZED SOUND PRESSURE SPECTRAL 
DENSITIES AT FREQUENCY F FOR THE THREE SPECTRAL SLOPES ARE: 


X1=CF/FC)**(SL1) 

X2=CF/FC)**CSL2) 

X3=CF/FC)**(SL3) 

DETERMINE THE ABSORPTION COEFFICIENT, ABCO IN DB/M, AT 
FREQUENCY F FOR AIR TEMPERATURE TK, RELATIVE HUMIDITY, 
HR, AND AIR PRESSURE OF 1.0 STD. ATM. 


ABCO=8. 686*CF#*2.)*CCC1. 84E-11)*C1. /PAD*SQRTCTK/TO)) 


1 +CCTK/T0)*#*®(-5./2.))*#CC1. 278E-2% 
2 CEXPC-2239.1/TK))/CFRO2+CCF#*2.)/FRO2))) 
3 +€0. 1068 *CEXPC-3352./TK))/CFRN2+CCF#%2.)/FRN2))))) 


THE ABSORPTION-LOSS FACTOR, ALF, OVER THE 300-M DISTANCE, 
IN THE DIRECTION OF PROPAGATION, IS: 


ALF=10.0**(-CABCO*300.0)/10.0) 


THE ARGUMENTS FOR THE NUMERATOR INTEGRAL OF THE 
BAND ATTENUATION FACTOR FOR THE THREE SPECTRAL SLOPES ARE: 


Y1CJUd )=X1 ALF 
Y¥2CJJ )=X2*ALF 
Y¥3CJJU)=X3 ALF 


CALL THE NUMERICAL INTEGRATION SUBROUTINE, QSF, TO INTEGRATE 
THE ARRAYS OF INTEGRAND VALUES OVER THE FREQUENCY RANGE OF THE 
IDEAL-FILTER PASSBANDS. 


CALL QSFCDELF,Y1,Z1,NN) 
CALL QSFCDELF,Y2,Z2,NN) 
CALL QSFCDELF,Y3,2Z3,NN) 


FIND BAND ATTENUATION FACTORS FOR THE THREE SLOPES. 
SLIPCK)=-10. 0*L0G10CCZ1(NN))/DEN1) 


SLONCK)=-10. 0*L0G10CCZ2(NN))/DEN2) 
SL12NCK)=-10.0*L0G10CCZ3(NN))/DEN3) 


260.000 C CE) CALCULATE ABSORPTION COEFFICIENTS BY SAE ARP866A 
261.000 C 

262.000 Bal. 328924-3.179768E-02*TC4+2.173716E-048TCHE2. 
263.000 1-1. 7496E-06*TC##3, 

264.000 ABSHUM=10. 0**(LOG10CHR)-B) 

265.000 AMOLMX=10. 0**(LOG10CFREQ(U))+8.42994E-03*TC-2. 755624) 
266.000 HMOLMX=SQRTCFREQCU)/1010.) 

267.000 HUMRAT =ABSHUM/ HMOLMX 

268.000 L=2 

269.000 IF CHUMRAT.LE.TABLEC(1))GOTO 100 

270.000 DO 70 M=3,57,2 

271.000 L=M+1 

272.000 1F CTABLE(M)-HUMRAT)70, 100, 80 

273.000 CONTINUE 

274.000 GOTO 100 

275.000 M=M~2 

276.000 IFCM.GE.3)GOTO 90 

277.000 M=M+2 

278.000 XA 1=HUMRAT-TABLE(M) 

279.000 XA O=HUMRAT-TABLE(M-2) 

280.000 XA 2=HUMRAT -TABLE (M42 ) 

281.000 X01=TABLE (M-2 )- TABLE (M) 

282.000 X02=TABLE (M-2 )-TABLE (M+2 ) 

283.000 X12=TABLECM)-TABLE (CM+2 ) 

284.000 ALPRAT=TABLECM-1 )®CXA1/X01)*CXA2/X02) 

285.000 -TABLE CM+1 )®CXA0/X01)*CXA2/X12) 

286.000 2 +TABLE (M43 )*(XA0/X02)*CXA1/X12) 

287.000 GOTO 110 

288.000 ALPRAT=TABLE(L) 

289.000 ALPMOL =ALPRAT *AMOLMX 

290.000 ALPCLA=10. 0**(2.05*LOGLOCFREQ(U)/1000.) + 1.1394E-03*TC 
291.000 1 - 1.916984) 

292.000 ABSORB=0.01*CALPCLA + ALPMOL) 

293.000 

294.000 THE ATTENUATION OVER A 300-M DISTANCE IS: 

295.000 

296.000 SAE 866(K )=ABSORB*300.0 

297.000 

298.000 CONTINUE 

299.000 

300.000 WRITEC6,130) FRQTICU), IFREQCU),CS126(K), SLIPCK), 
301.000 1ISL6NCK), SLI2NCK), SAE 866(K),K=1,5) 

302.000 FORMAT(T1,A1,16, 5(5F5.1)) 

303.000 C 

304.000 CONTINUE 

305.000 

306.000 WRITEC6, LSODFRQTIC25) 

307.000 FORMAT(T1,A1) 

308.000 WRITEC6, 160) FRQTIC26) 

309.000 FORMAT(T1,A1,T17,'60 PCT RH',16X,'70 PCT RH! 
310.000 1'80 PCT RH',16X,'90 PCT RH',15X,'100 PCT RH'§ 
311.000 WRITEC6,170) FRQTIC27) 

312.000 FORMATCT1,A1,78,5' CA) CB) CC) (CD) CED") 
313.000 C 


16X, 


a aa weenie 


314.000 
315.000 
316.000 
317.000 
318.000 
319.000 
320.000 
321.000 
322.000 
323.000 
324.000 
325.000 
326.000 
327.000 
328.000 
329.000 
330.000 
331.000 
332.000 
333.000 
334.000 
335.000 
336.000 
337.000 
338.000 
339.000 
340.000 
341.000 
342.000 
343.000 
344.000 
345.000 
346.000 
347.000 
348.000 
349.000 
350.000 
351.000 
352.000 
353.000 
354.000 
355.000 
356.000 
357.000 
358.000 
359.000 
360.000 
361.000 
362.000 
363.000 
364.000 
365.000 
366.000 
367.000 
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DO 260 J=1, 24 
HR=50.0 

DO 240 K=6,10 
HR=HR+10.0 


CALCULATE THE EXACT BAND CENTER FREQUENCIES FROM THE 
INTERNATIONAL STANDARD BAND NUMBERS, ISBN. 


ISBN=J+16 
FC=10.0**(1SBN/10.0) 


THE EXACT LOWER, FL, AND UPPER, FU, BANDEDGE FREQUENCIES FOR 
THE FILTER BAND ARE: 


FL=FC/SQRTCRF) 
FU=FC*SQRTCRF) 


THE WIOTH OF THE EQUALLY SPACED FREQUENCY INTERVALS, DELF, 
OVER THE PASSBAND OF THE FILTER IS: 


DELF =(FU-FL)/STEPCuU) 


THE VALUE OF THE FREQUENCY AT ONE STEP BELOW THE LOWER 
BANDEDGE FREQUENCY IS: 


F=FL-DELF 


THE TOTAL NUMBER OF FREQUENCY STEPS IN THE NUMERICAL 
INTEGRATION OVER THE PASSBAND OF THE IDEAL FILTER IS: 


NN=1+STEPCuU) 

CALCULATE THE DENOMINATOR TERMS USED IN THE INTEGRATION 
METHOD TO DETERMINE THE BAND ATTENUATION OVER EACH OF THE 
THREE SPECTRAL SLOPES. 

THE DENOMINATOR TERM FOR THE +1.0 DB/BAND SLOPE IS: 


DENI=FC#(1.0/CSL141.0))*CCRF**#(CSL141.0)/2.0))- 
1 CRF*#*(-(SL141.0)/2.0))) 


THE DENOMINATOR TERM FOR THE -6.0 DB/BAND SLOPE IS: 


DEN 2=FC#(1.0/CSL 241.03) *CCRF¥*CCSL 241 .09/2.0))- 
1 CRF#*%(-(SL241.0)/2.0))) 


THE DENOMINATOR TERM FOR THE -12.0 DB/BAND SLOPE IS: 


DENS=FC#C1. 0/CSL341.0))*CCRF*#*CCSL341.0)/2.0))- 
1 CRF##(-(CSL34+1.0)/2.0))) 
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369.000 
370.000 
371.000 
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CA) CALCULATE ABSORPTION COEFFICIENTS BY S1.26-1978 


LPSOP0=10.79586*(1.-(TO1/TK)) - 5.02808*LOG1OCTK/TO01) 


1 +1. 504746 ~42¢(1.-10.%*(-8. 29692*CCTK/TO1)-1.))) 
2 +0.42873E~3*C10.%*C4.76955%C1.-CTOL/TK)))-1.) 
3 -2.2195983 


PSOP0=10.**LPSOPO 

H=HR*PSOP0/PA 

FRO2=PA¥( 24.44. G16 04*H#®CC0.054H)/C0.391+H))) 
FRN2=CPA/SQRTCTK/T0))*(9.+ 


1 350.*%H*EXP(-6.142*C(CCTK/TODR*®(C-1./3.))-1.)9)) 
ALPHA=CFC#*2.)*CCC1. 84E-11)%(1./PAD*®SQRTCTK/TO)) 

1 +CCTK/TOD#*(-5./2.))*#CC1. 278E-2% 

2 © CEXPC=-2239.1/TK))/CFRO2Z+CCFC#*2.)/FROZ))) 

3 +€0.1068*CEXP(-3352./TK))/CFRN24(CFC#*2.)/FRN2))))) 


THE ABSORPTION COEFFICIENT, AFC IN DB/M, AT THE EXACT 8AND 
CENTER FREQUENCY BY S1.26-1978 IS: 


AFC=8.686*ALPHA 
THE ATTENUATION OVER A 300-M DISTANCE IS: 
$126CK)=AFC*300.0 


CB), €C), AND (CD) CALCULATE BAND ATTENUATION FACTORS BY 
THE INTEGRATION METHOD FOR THE THREE SPECTRAL SLOPES. 


DO 180 JJ=1,NN 


DETERMINE THE FREQUENCIES, F, AT EACH STEP OVER THE 
FREQUENCY RANGE OF THE FILTER PASSBAND. 


F=F+DELF 


THE VALUES OF THE NORMALIZED SOUND PRESSURE SPECTRAL 
DENSITIES AT FREQUENCY F FOR THE THREE SPECTRAL SLOPES ARE: 


X1LaCF/FC)*#¥CSLI) 

X2=CF/FC)**CSL2) 

X3=CF/FC)##CSL3) 

DETERMINE THE ABSORPTION COEFFICIENT, ABCO IN DB/M, AT 
FREQUENCY F FOR AIR TEMPERATURE TK, RELATIVE HUMIDITY, 
HR, AND AIR PRESSURE OF 1.0 STD. ATM. 


ABCO=8. 686 *CF#%2.)#CCC1. 84E-11)*C1. /PA*®SQRTCTK/TO)) 


1 +CCTK/TOD®#C25./2.))*CC1. 278E-2% 
2 CEXP(=-2239. 1/ TKI /CFRO24+CCFR*2.)/FROZ))) 
3 +#C€0. 1068*CEXPC-3352./TKII/CFRN2+CCFR82.)/FRN2))))) 


ee eee ee 


420.000 
421.000 
422.000 
423.000 
424.000 
425.000 
426.000 
427.000 


“428.000 


429.000 
430.000 
431.000 
432.000 
433.000 
434.000 
435.000 
436.000 
437.000 
438.000 
439.000 
440.000 
441.000 
442.000 
443.000 
444,000 
445.000 
446.000 
447.000 
448.000 
449.000 
450.000 
451.000 
452.000 
453.000 
454.000 
455.000 
456.000 
457.000 
458.000 
459.000 
460.000 
461.000 
462.000 
463.000 
464.000 
465.000 
466.000 
467.000 
468.000 
469.000 
470.000 
471.000 
472.000 
473.000 


olekene) aan 


(omens) AAAAN 


fomeke) 


THE ABSORPTION-LOSS FACTOR, ALF, OVER THE 300-M DISTANCE, 


IN THE DIRECTION OF PROPAGATION, IS: 


ALF =10.0**(-CABCO*300.0)/10.0) 
THE ARGUMENTS FOR THE NUMERATOR INTEGRAL OF THE 


BAND ATTENUATION FACTOR FOR THE THREE SPECTRAL SLOPES ARE: 


1806 


Y1CJuU)=X1 ALF 
Y2CJU )=X2 ALF 
Y3CJUJU)=X3 “ALF 


CALL THE NUMERICAL INTEGRATION SUBROUTINE, QSF, TO INTEGRATE 


THE ARRAYS OF INTEGRAND VALUES OVER THE FREQUENCY RANGE OF THE 
IDEAL-FILTER PASSBANDS. 


190 
200 


CALL QSFCOELF,Y1,Z1,NN) 
CALL QSFCOELF,Y2,Z2,NN) 
CALL QSFCDELF,Y3,Z3,NN) 


FIND BAND ATTENUATION FACTORS FOR THE THREE SLOPES. 


SLIP(K)=-10.0*L0G10CCZ1(NN))/DEN1) 
SLONCK )=-10.0*L0G10CCZ2(NN))/DEN2) 
SL12NCK )=-10.0*L0G10C(Z3(NN))/DEN3) 


CE) CALCULATE ABSORPTION COEFFICIENTS BY SAE ARP866A 


B=1.328924-3.179768E-02*TC+2.173716E-O4*TC#EQ. 
1-1.7496E-06*TC##3, 

ABSHUM=10.0**(LOG10(HR)-8) 
AMOLMX=10.0**CLOG10CFREQCU))+8.42994E-03*TC~2.755624) 
HMOLMX=SQRTCFREQ(CU)/1010.) 

HUMRAT =ABSHUM/ HMOLMX 

L=2 

IF CHUMRAT.LE.TABLEC1))GOTO 220 

DO 190 M=3,57,2 

L=M+1 

IF CTABLECM)=HUMRAT 190, 220, 200 

CONTINUE 

GOTO 220 

M=M-2 

IFCM.GE.3)GOTO 210 

M=M+4+2 

XA L=HUMRAT-TABLE(M) 

XA 0=SHUMRAT-TABLE (M-2 ) 

XA 2=HUMRAT -TABLE (M42 ) 

XO01=TABLE(M-2 )-TABLECM) 
X02=2TABLE(M-2 )=-TABLE (M42 ) 
X12=TABLECM)-TABLE(M42 ) 
ALPRAT=TABLE CM-1 0#®(XA1/X01)#CXA2/X02) 

1 -TABLECM+1 )*(XA0/X01)*CXA2/X12) 

+TABLE CM+3 )®(XA0/X02)*CXA1/X12) 

GOTO 230 
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ALPRAT=TABLEC(L) 

ALPMOL =ALPRAT ®AMOLMX 

ALPCLA=10. 0**(2.05*LOGIOCFREQ(U)/1000.) + 1.1394E-03*TC 
1 - 1.916984) 

ABSORB=0.01*CALPCLA + ALPMOL) 


THE ATTENUATION OVER A 300-M DISTANCE IS: 
SAE 866 CK J=ABSORB*300.0 
CONTINUE 


JJ=JU+27 

WRITEC6, 250) FROQOTICUu), IFREQCU),CS126(K), SLIPCK), 
ISLONCK),SLI2NCK), SAE866(K),K=6, 10) 
FORMAT(T1,A1, 16, 5(5F5.1)) 


CONTINUE 
CONTINUE 
END 


SUBROUTINE QSFCH,Y,Z,NDIM) 


aeee 
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SEIS BSS ss 8288686 36 38 38 28 38 58 38 58 28 38 SE 38 BE TC 
SUBROUTINE QSF FROM I8M SCIENTIFIC SUBROUTINE PACKAGE (SSP), 
FIFTH EDITION,AUGUST 1970. 


PURPOSE 
TO COMPUTE THE ARRAY OF INTEGRAL VALUES FOR A GIVEN 
EQUIDISTANT TABLE OF FUNCTION VALUES. 


USAGE 
CALL QSF CH,Y,Z,NDIM) 


DESCRIPTION OF PARAMETERS 
H -THE INCREMENT OF ARGUMENT VALUES. 
Y -THE INPUT ARRAY OF FUNCTION VALUES. 
Zz -THE RESULTING ARRAY OF INTEGRAL VALUES. Z MAY BE 
IDENTICAL WITH Y. 
NDIM <THE DIMENSION OF ARRAYS Y AND Z. 


REMARKS 
NO ACTION IN CASE NDIM LESS THAN 3. 
NO SUBROUTINES OR FUNCTION SUBPROGRAMS ARE REQUIRED. 


METHOD 
BEGINNING WITH Z(1)=0, EVALUATION OF ARRAY Z IS DONE BY 
MEANS OF SIMPSON'S RULE TOGETHER WITH NEWTON'S 3/8 RULE OR 
A COMBINATION OF THESE TWO RULES. TRUNCATION ERROR IS OF 
ORDER H¥#5 C1.E.,FOURTH ORDER METHOD). ONLY IN CASE NOIM=3 
IS THE TRUNCATION ERROR OF 2(2) OF ORDER H*#4, 


38 28 56 26 36 56 56 28 BE OE 9 28 26 26 SE 28 36 36 9E BE OE 26 36 36 BE SE 28 98 BE 98 98 OE 96 38 E38 25 98 BE OE SE 38 98 38 HE 38 OE OE OE SE SE OE OE BE aE Oe oe Se See 
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530.000 C 

531.000 DIMENSION Y(1),2(€1) 

532.000 C 

533.000 HT=0.3333333H 

534.000 IFCNDIM~5)7,8,1 

535.000 C 

536.000 C NOIM IS GREATER THAN 5. PREPARATION OF INTEGRATION LOOP. 
537.000 C 

538.000 1 SUM1=Y(2)+YC2) 

539.000 SUM1=2SUM1+SUM1 

540.000 SUMISHT*CYC1)+SUM1+Y(3)) 
541.000 AUX1=YC4)+Y C4) 

542.000 AUX1=AUX1+AUX1 

543.000 AUX1=SUM14HT *CY C3 )4AUX1+Y(5)) 
544.000 AUX2=HT#CY C143 .8752CY C2947 C5 42.625 #CV C35 +¥C4))4+Y(6)) 
545.000 SUM2=Y(5)+Y(5) 

546.000 SUM2=SUM2+SUM2 

547.000 SUM 2=AUX2-HT *CY C4 +SUM2+Y (6 )) 
548.000 Z€1)=0.0 

549.000 AUX=Y(3)+Y(3) 

550.000 AUX =AUX+AUX 

551.000 ZC 2 )=SUM2-HT *CY C2 )+AUX4Y(4)) 
552.000 Z(3)=SUM1 

553.000 ZC4)=SUM2 

554.000 IFCNDIM-6)5, 5,2 

555.000 

556.000 INTEGRATION LOOP. 

557.000 

558.000 DO 4& I=7,NDIM, 2 

559.000 SUM1=AUX1 

560.000 SUM2=AUX2 

561.000 AUX1=YC1I-1)4YCI-1) 

562.000 AUX1=AUX14AUX1 

563.000 AUX1=2SUM14HT#CY CI -2 4AUX14YCI)) 
564.000 ZC 1=-2)=SUM1 

565.000 IFCI-NDIM)3,6,6 

566.000 AUX2=YC1)+YC1) 

567.000 AUX 2=AUX2+4AUX2 

568.000 AUX2=SUM24HT CY CI ~1 )+4AUX24+Y C141 )) 
569.000 ZCI -1)=SUM2 

570.000 ZCNDIM-1 )=AUX1 

571.000 ZCNDIM)=AUX2 

572.000 RETURN 

573.000 ZCNDIM-1 )=SUM2 

574.000 ZCNDIM)=AUX1 

575.000 RETURN 

576.000 

577.000 END OF INTEGRATION LOOP FOR NDIM GREATER THAN 5. 
578.000 

579.000 7 IF CNDIM-3)12,11,8 

§80.000 C 

581.000 C NOIM IS EQUAL TO & OR 5. 
582.000 C 
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| | 583.000 8 = SUM2=1.125*HT#CYC1)+¥(2 +72 4¥ C241 C3 41 C3 4¥(3 49 C4)) 


584.000 SUM1=Y(2)+Y(2) 
585.000 SUM 1=5SUM1+SUM1 
586.000 SUM 1SHT*CY (1 )+SUM1+Y(3)) 
587.000 2€1)=0.0 
588.000 AUX1=Y(3)+Y(3) 
589.000 AUX1=AUX1+AUX1 
§90.000 2€2=SUM2=HT ¥CY C2 +AUX1+Y(4)) 
591.000 IFCNDIM-5)10,9,9 
592.000 C 
593.000 C HERE NDIM IS EQUAL TO 5. ; 
§94.000 C 
$95.000 9 AUX1=Y(4)+Y (4) 4 
596.000 AUX1=AUX1+4AUX1 ; 
597.000 ZC5 D=SUM14HT ®CY¥ C3 )4+AUX14+YC5 )) i 
598.000 C 
599.000 C HERE NDIM IS EQUAL TO 4. i 
600.000 C ; 
601.000 10 2Z¢€3)=SUM1 
602.000 Z(4)=SUM2 i 
603.000 RETURN i 
604.000 C i 
605.000 ¢ HERE NDIM IS EQUAL TO 3. 
606.000 C t 
607.000 11 SUMI=SHT#C1L, 25*YC1)+Y(2)+Y(2)-0. 25%Y(3)) 
608.000 SUM2=¥(2)+Y(2) 
609.000 SUM 2=SUM2+SUM2 
610.000 Z€3 J=HT*®CY C1 )+SUM2+Y (3 )) j 
611.000 2(1)=0.0 
612.000 Z(2)=SUM1 
613.000 12 RETURN 
614.000 END ' 
| 
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